Fluctuating Filaments Under Tension - From Flexible Chains to Rigid Rods 
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We develop a novel biased Monte-Carlo simulation technique to measure the force-extension curves 
and the distribution function of the extension of fluctuating filaments stretched by external force. 
C^) ' The method is applicable for arbitrary ratio of the persistence length to the contour length and for 

arbitrary forces. The simulation results agree with analytic expressions for the force-extension curves 
and for the renormalized length-scale-dependent elastic moduli, derived in the rigid rod and in the 
strong force limits. We find that orientational fluctuations and wall effects produce non-Gaussian 
| distributions for nearly rigid filaments in the small to intermediate force regime. We compare our 

. results to the predictions of previous investigators and propose new experiments on nearly rigid rods 

Pjr^ ■ such as actin filaments. 

VO ' PACS numbers: 87.15.Aa,87.15.Ya,05.40.-a 
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I. INTRODUCTION 

o 

^ . Recent experimental studies of the mechanical properties of biological materials, ranging from individual DNA 
molecules^ to single actin filaments and their networks 00. El- led to revival of interest in the worm-like chain model 
of semiflexible biopolymers. In particular, the force-extension curves predicted by this model in the long chain limit 
(i.e., when the length of the chain is much larger than the persistence length) were shown to agree quantitatively with 
| experimental results on dsDNA at large elongations p|. However, many interesting questions concerning worm-like 
c"j . chains under tension are still open and there are delicate issues concerning even the linear response of the worm-like 
chain to stretching forced . Also, while the equilibrium distribution function of the extension was calculated for the 
worm-like chain|J0, Q and fluctuating ribbon0] models, only approximate (u sing the Gaussian approximation) 
calculations of the distribution function for stretched chains were reported to date|ll|. Such distributions have been 
7— ( | measured in experiments on stretched DNA molecules anc l provide a much more stringent test of the theory than 
^ ■ the force-extension curves. Finally, even though all mechanical experiments involve biopolymers attached to surfaces, 
none of the present theories incorporates wall effects on force-extension curves and on the corresponding distributions. 
While such effects are negligible for nearly fully stretched chains, they are expected to be important in the weak force 
; regime. 

In previous work |l(A Il3j . we have discussed how to generate via Monte-Carlo equilibrium ensembles of fluctuating 
• ribbons with spontaneous curvature and/or twist. This method was based on using the intrinsic geometry of the 
filament, in which representation the energy is purely quadratic. This enables a direct generation of independent 
samples distributed according to the Gibbs equilibrium ensemble. When external forces act on the filaments, the 
energy is no longer a quadratic function of the curvatures, and the problem is now not only nonlinear but also 
nonlocal. In this paper we discuss a biasing method which succeeds in efficiently generating Monte-Carlo samples for 
the simplest case of an external force, namely that of the stretched filament. This problem was studied by Marko 
and Siggia, (MS), 0,01 an d others who calculated the force-extension curve for long worm-like chains. We 

will reproduce these results as well as extend the analysis to the case of nearly stiff filaments with length comparable 
£j , to or shorter than the persistence length, and in addition calculate the distribution of extensions for various forces. 
Lastly, we will show how we can also easily incorporate steric constraints; we will demonstrate this in the context of 
a restriction of the filament to a half-space by the presence of a wall. 

In sections 2 and 3 we review the general theory of fluctuating ribbons (we define a ribbon as a filament with an 
asymmetric cross-section) and the Frenet algorithm for computer simulations of such ribbons, respectively. In section 
4 we describe the biased Monte-Carlo simulations of filaments under tension and present the resulting force-extension 
curves and distribution functions. Wall effects are considered in section 5. In section 6 we discuss our results and 
outline directions for future research. In Appendix A we expand the Euler angles in terms of the principal curvatures 
and then use this expansion in Appendix B in order to calculate the partition function of a stretched filament. In part 
1 of this appendix we derive an analytic expression for the partition function and obtain a force-extension relation for 
arbitrary force in the limit of large persistence length (expansion about the rigid rod limit). In part 2 of the appendix 
we derive asymptotic expressions for the above quantities in the limit of large force (for arbitrary rigidity). Finally, 
in Appendix 3 we present the calculation of the renormalization of the elastic moduli of the various Fourier modes of 
the curvature. 
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II. GENERAL THEORY 

The linear elasticity of ribbons is described by the elastic energy p"7] 

1 3 , r L 

E e i = r 22 bk / ds ( Uk ~ Wofe ) 2 ' W 
1 fe=i 70 

where s is the distance measured along the contour of the ribbon, the coefficients b\ and 62 are the moduli associated 
with bending along the two principal axes of inertia of the cross section (they differ if the cross section is not 
circular), and 63 is the twist modulus. In this paper we treat {bi} as given material parameters of the ribbon. The 
functions {cuk(s)} and {wot(s)} are the generalized curvatures in the deformed and the stress-free states of the ribbon, 
respectively. These curvatures completely determine the three dimensional conformation of the ribbon through the 
generalized Frenet equations 

ds = - 2^£ijkUj{s)t k (s). (2) 

Here t3 is the unit tangent to the centerline and the unit vectors ti and t2 are oriented along the principal axes of 
the cross section (sijk is the antisymmetric unit tensor). Note that since these equations describe a pure rotation of 
the triad {t^s)} as one moves along the contour of the ribbon, the constraint |t 3 | = \dr(s)/ds\ = 1 is automatically 
satisfied in this intrinsic coordinate description. 

The physical frame of the ribbon {tj(s)} can be related to the original Frenet description of space curves in terms 
of a unit tangent (which coincides with t 3 ), normal n and binormal b, provided that in addition to the Frenet triad 
{b(s), n(s), ts(s)} associated with the centerline of the ribbon, one introduces an angle a(s) that describes the rotation 
of two of the principal axes of the cross section t\, t2 with respect to the binormal b and normal t (see Fig. 1 in the 
second of Re fs. 0] ), 

ti = b cos a + n sina, t2 = — bsina + ncosa. (3) 
The three dimensional configuration of the centerline is determined by the Frenet equations 

db dn dt 3 

— = rn, — = Kt a - rb, — = -sn (4) 

ds ds ds 

where the curvature k(s) and torsion t(s) are treated as known functions of position along the contour. The three 
generalized curvatures {^(s)} , can be expressed in terms of these parameters and the angle a(s), 

(jJi = kcosq!, u>2 = Ksina, u>3 = r + da/ds. (5) 



III. THE FRENET ALGORITHM 



Since the energy is a quadratic form in the deviations 5u>k of the generalized curvatures from their values in the 
stress free state, Eq. itlj) is valid as long as the characteristic length scale of the deformation is much larger than the 
diameter of the ribbon[18|. When a ribbon undergoes fluctuations in the presence of a thermal bath at temperature 
T, this energy determines the statistical weight of the configuration {uJk}- The statistical average of any functional of 
the configuration A{{uik\) is defined as the functional integral [l7| 

{A ({W * })> = jDiS^e-^^T • (6) 

where k is the Boltzmann constant. From the form of the energy, Eq. ^ it follows that the distribution of Suik is 
Gaussian, with zero mean ((<5u>j(s)) = 0) and second moments given by 



(SuiWSwjtf)) = — S ZJ 5(s - s'). (7) 

In order to study the distribution functions of the various fluctuating quantities (e.g., the extension i?), we developed 
what we entitled the Frenet algorithm, which is based on the combination of statistical methods and differential 
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geometry and is described in the following. Examination of Eq. (JJJ shows that the decoupling of the "noises" {Suj^} 
in the intrinsic coordinate representation permits an efficient numerical generation of independent samples drawn 
from the exact canonical distribution. The Gaussian distribution of the Suik means that each 5(jJk(s) can be directly 
generated as one of a string of independent random numbers drawn from a distribution symmetric about the origin 
with width •v/rcT/(6fcAs) where As is the discretization step length. The remaining task is to construct the curve 
using the Frenet equations with Wfc(s) = oJok{s) + 5uJk{s). The Frenet equations are best integrated by stepping the 
basic triad tk forward in s through a suitable small rotation. In this way, the orthonormality of the triad is guaranteed 
to be preserved up to machine accuracy. To construct this rotation matrix, we begin with Eqs. J2J and, defining the 
three vectors v x = (tf,t^,tf), and so forth, we can write this equation as 

(8) 

As 

where A is an antisymmetric matrix with elements Aij = J^ fc EijkU>k- We now discretize Eq. (jSJ) as 

i/(s + As) = Oi/(s), (9) 



0=(1 + ^a)(i-^A) . (10) 



where the orthogonal matrix O is 

n = (\ + : 

2 J \ 2 

The choice of discretization step length (L/As is the number of points) depends on the amplitude of the random noise, 
y/kT/(bkAs). We found that 400 points gave sufficient accuracy (i.e., was limited by the statistical noise determined 
by the number of independent samples) for the calculation of the distribution function of the extension, when the 
smallest of the moduli satisfied kTL/bf in < 1. 

This algorithm is optimal for studying the fluctuations of a free ribbon. However, if the ribbon is acted upon by 
external forces and/or constraints, the algorithm must be modified. This is described in the next section. 

IV. BIASING THE MONTE-CARLO 

Consider a filament with a circular cross-section (b\ = 62) which is characterized by a bending persistence length 

The simplest example of a filament under tension is that of a constant force acting on one end and directed along 
the x-axis, with the second end constrained to lie at the origin but otherwise free to rotate about it. The energy of 
the filament thus picks up an additional contribution —Fx, where x is the projection of the end to end vector on the 
direction of the force. The problem is that a; is a complicated nonlinear function of all the Su>, and so the external 
force destroys the Gaussian structure of the energy on which the generation of the Su> depends. One can simply 
incorporate the extra contribution to the energy by treating it as a weighting factor, or bias, and suitably normalizing 
all averages. Thus, for example, one can calculate the mean value of x via 

^=5^ (12) 

where the sum is over configurations generated via the Frenet algorithm described above and / = F/kT is the force in 
units of inverse length. For small forces, this works well. However, as the force increases, the probability of generating 
highly extended configurations which are heavily favored by the bias, is so small that the sampling becomes ineffective. 
One therefore needs a way to bias the generation of configurations towards those that are extended, so they can be 
sampled appropriately. 

This task is not simple, since the connection between the Jw's through which the filament is parameterized and 
the value of x is very complicated. We can gain insight into this by transforming to Fourier space and looking 
at the expectation value of the moments of the n'th Fourier component of the curvature, Suj(n), using the above 
algorithm, Eq. ^1 for small to intermediate forces where it is still effective. For an extensional force, we do not 
expect spontaneous curvature to develop, and indeed the first moments of the Fourier components all vanish. The 
second moments, however, do exhibit an interesting structure, as seen in Fig. ^ where the expectation value of 
the squared amplitude of the low n Fourier modes are presented for the case a — 3, L — 2tt. For F = 0, the 
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second moments are independent of n, as expected. When the filament is subjected to a tensile force we see that 
the small n modes are suppressed, with the degree of suppression increasing with F and decreasing with n. In the 
case of large a, where the fluctuations are in any case small, we observe that only the n = mode is significantly 
affected by the force (not shown). The suppression of the various curvature modes can be calculated analytically 
when fluctuations are small, namely either in the large FL limit, or for any F provided \ = \J FL 2 /(AakT) <C 1 
(see Appendix). The calculation confirms the above-mentioned trends. We learn from this that the high Fourier 
modes (or equivalently the small-scale fluctuations) of the filament are essentially unaffected by the force, whereas 
the long-wavelength fluctuations are significantly damped. 

This suggests the possibility of biasing the Monte-Carlo by selectively dampening the small-n modes. The idea 
is to generate the <5u?s in Fourier space, with appropriate amplitudes for the various modes, and construct the real- 
space 8uj(s) via a Fast Fourier Transform. Of course, this bias has to be corrected explicitly in the weighting in the 
ensemble average. We do this in a self-consistent fashion, choosing the amplitude of each Fourier mode to reproduce 
the measured square magnitude of the mode. In practice, it is sufficient to do this only for the first few modes, which 
are the ones most severely impacted by the force. Also, in order to reduce the noise in the measurement of the mode 
amplitudes, we only calculate the averages after every block of N/1Q realizations and update accordingly, where N 
is the total number of realizations we intend to generate for these parameters. The force dependence of the mode 
amplitudes is obtained by a stepwise increase of F, where for each value of the force we use as a starting point the 
amplitudes measured at the previous value. 

There is one other degree of freedom of the problem that also must be treated appropriately, namely the initial 
orientation of the filament at its constrained end. This is because the force not only straightens the filament, but 
also orients it toward its direction. We have employed two different algorithms to deal with this problem. The 
simpler approach is to simply bias the choice of initial orientation, and correct for this in the weighting. The zero 
force distribution of the orientation angle, 9, of the tangent vector with respect to the direction of the force is 
P(9) = sin9d9 = d(cos6). In the presence of the force we modified this to P(cos9) — exp(a cos 9) d(cos 9), where one 
can choose a to reproduce the measured mean of cosQ. 

A more sophisticated approach is to explicitly perform the integral over 9. We generate a given configuration and 
calculate its end-to-end distance R and angle 9 cn d that the tangent to the end of the filament makes with the direction 
of the end-to-end vector R. In order to evaluate the contribution of this realization to the average elongation, Eq. 
1121 one can perform a rigid rotation of this configuration about the end point and calculate the contribution to the 
statistical sum, from each angle of rotation 9. The resulting integrals over 9 can be expressed in terms of R and 9 cn d 
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FIG. 2: Extension along the force direction of a filament, as measured by our Monte-Carlo procedure, together with the 
Marko-Siggia interpolation formula and our analytic approximation, Eq. (IB23t for L — 2ir. Left: a = 0.75; Right: a = 48. 
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where the curly braces indicate the integrand for the numerator and denominator of the ensemble sum, respectively. 
These integrals can easily be evaluated numerically, using standard techniques. This second method is more efficient, 
but more involved to implement. 

Using these techniques, we have measured the force-extension curve x vs. /. We present the results for a = .75, 
L = 2ir in the left box of Fig. [5] together with the interpolation formula of MS 



f " = I + 4(1 -\,LY I ,l4! 



We see that for this value of a (a <C L) the agreement is very good. In the right half of Fig. [21 we present the results 
for a = 48, again with L = 2ir. Here the MS formula does not work for small and intermediate values of the force, 
since we are not in the long-chain limit (L ^> a) for which it was derived. On the other hand, for this case analysis can 
also be done, as presented in the Appendix, and the analytic formula derived there is in perfect agreement with the 
simulations (see the right box of Fig. Notice that the MS expression systematically overestimates the extension 
of the filament and approaches our results only in the infinite force limit. The physical basis of the deviation is the 
large orientational entropy which is not accounted for in the MS calculation. 

Inspection of the right box of Fig. [3 suggests that both the simulation and our analytical approximation (Eq. IB 121) 
yield a finite initial slope of the extension versus force curve in the rigid rod limit (a» L). Since the effective spring 
constant is proportional to the inverse of the initial slope of this curve this suggests that the linear response of the 
rigid rod to small stretching forces does not depend on the persistence length a and is a function of its length L only. 
The above expectation is confirmed by calculations reported in Appendix B where we show that the corresponding 
spring constant is 3kT/L 2 . Obviously, this behavior differs from the flexible chain limit (a <C L) for which the spring 
constant scales as^ kT/(aL). It also differs from the predictions of Ref. 0] where it was argued that the effective 
spring constant diverges in the rigid rod limit, as akT/L 3 . Notice that such a divergence is not expected for a single 
filament under tension since there is a finite contribution to the susceptibility (and thus to the response to the external 
force) due to orientational fluctuations, even in the rigid rod limit. The difference between our result and that of Ref. 
] stems from the fact that the latter considers a network in which filaments are not free to rotate and the only 
angular averaging is over the randomly distributed orientations of different filaments in the network. 

Not only can we calculate the average x, but we can also measure the entire distribution P(x). We show this in Fig. 
13 for a — 0.75 and 6. In the former case, the persistence length is much shorter than length of the filament and the 
force-free distribution is nearly Gaussian. The peak of the distribution moves to higher values of the extension and its 
width decreases with increasing force. In the latter case (persistence length approximately equal to the length of the 
filament), the force-free distribution is nearly flat up to x ~ ±5, as expected for random sampling of orientations by 
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FIG. 3: Distribution of the extension along the force direction, as measured by our Monte-Carlo procedure, for the cases / = 0, 
1, 2 and 3, with L = 2ir. Left: a = 0.75; Right: a — 6. 

a nearly rigid rod (due to thermal fluctuations, the effective length of this rod is somewhat shorter than its contour 
length, 2ir). The application of small force (in the range F/kT < 1) orients the filament and leads to dramatic change 
of the distribution. At larger forces, the distribution changes slowly with the force, in a manner similar to the case 
a = 0.75. The transition from orientation-dominated to stretching-dominated behavior of nearly rigid rods is clearly 
observed in the right box of Fig. |2J where pronounced deviations from the Marko-Siggia formula are observed in the 
weak force regime (orientational effects are not taken into account in the MS theory). 

V. STERIC CONSTRAINTS 

It is fairly straightforward to incorporate steric constraints into our calculations. As an example, we consider the 
case where the filament is confined to the upper half space x > 0. After generating each configuration according to 
the procedure outlined in the previous section, we simply reject any configuration which violates the constraint. Our 
method of integrating over the initial orientation angle 9 can also be employed. We have simply to compute the limits 
of the integral over 6 which can be done straightforwardly given the configuration. We present the results for the 
force-extension curve in Fig. 0]and f° r the distribution P(x) in Fig. [3] We see that, as one might expect, the effect of 
the constraint is most pronounced at small F, where the straightening/ aligning of the curve is smallest. For vanishing 
F, the average extension is positive as a result of the constraint. The constraint basically carves out a "hole" in the 
distribution at small x, resulting in a very asymmetric distribution for small force. For F/kT > 1 the distribution 
function is practically unaffected by the presence of the wall and coincides with that shown in Fig. [3J 

VI. DISCUSSION 

In this work we used biased Monte-Carlo simulations, based on an extension of our Frenet algorithm, to calculate 
the force-extension relations and the distribution functions of the extension for filaments of arbitrary persistence to 
contour length ratio. We showed that in the presence of external force, curvature fluctuations are no longer Gaussian 
and that their long wavelength Fourier modes are progressively suppressed by the force. This effect was taken into 
account in a self-consistent manner in the simulations by gradually increasing the force and adjusting the corresponding 
Fourier amplitudes of curvature fluctuations. Analytic results were obtained both for small fluctuations about the 
rigid rod limit, and in the limit of large force. In the former case, we calculated the force-induced renormalization 
of the length-scale-dependent bending moduli that govern the various Fourier modes of curvature. As expected, the 
largest effect is on the zeroth (constant curvature) mode; the corresponding modulus increases quadratically and 
linearly with the force in the weak and strong force limits, respectively. The underlying physical picture is reminiscent 
of the suppression of excluded volume by tensile force in the self-avoiding random walk model of polymers - while 
large scale fluctuations are strongly suppressed, small scale fluctuations (within the so called Pincus blobs |2(j) are 
not affected by the force. 
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FIG. 4: Force-extension curve with and without the presence of a wall confining the filament to x > for L — 2ir, measured via 
our Monte-Carlo procedure. Left: a = 0.75, together with the Marko-Siggia formula; Right: a = 6, together with our formula, 
Eq. (|B2SJ. 




FIG. 5: Distribution of the extension along the force direction with and without the constraining wall, for / = 0, 1, with 
L = 2-7T. Left: a — 0.75; Right: a = 6. Distributions at higher forces are essentially unaffected by the constraint. 



For flexible filaments with a <ti L our simulation reproduces the Marko-Siggia force-extension curves. However, 
our simulations and analytic results show that as the rigid rod limit a > L is approached, the MS expression 
increasingly overestimates the response of the filament at finite values of the force, and approaches the correct result 
only asymptotically, in the limit of infinite force. The discrepancy can be attributed to the neglect of orientational 
effects (the constant curvature mode is not taken into account) in the MS derivation which is strictly valid only 
in the a/L — > limit. These orientational effects are however seen in the short-chain limit of the Gaussian model 
investigated recently by Winkler 11]. The large force limit of the Gaussian model for short chains is however very 
different from the worm-like chain model studied herein. For long chains, the Gaussian model can reproduce cither 
the weak or the strong force limit, by appropriate choice of the persistence length, but not both. This is an outcome 
of the length-scale dependence of the effective moduli discussed above and in Appendix C. 

Examination of the distribution functions for the extension shows that increasing rigidity leads to a more rapid shift 
of the peak of the distribution and narrowing of its width, with external force. The largest deviations from Gaussian 
statistics of the distributions are observed for stiff filaments in the weak to intermediate force range (F/kT < 1), where 
fluctuations of initial orientation are not yet completely suppressed and where wall effects play an important role. 
While it may be difficult to approach the rigid rod limit in mechanical experiments on dsDNA that has a persistence 
length of 500 A, it can be readily accessed by experiments on actin filaments 2] for which the persistence length of 
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order 10 /im. 

Finally, we would like to comment on possible extensions of the present work. Since we only considered here fila- 
ments with symmetric cross section and without spontaneous curvature, twist rigidity did not affect the distribution 
of extension and played no role in our considerations. However, since twist rigidity affects even the equilibrium distri- 
butions of ribbons with strongly asymmetric cross sections|10|. it is expected that both twist rigidity and spontaneous 
twist will affect the deformation of such ribbons as well. Perhaps even more interesting is the effect of spontaneous 
curvature. Indeed, we have shown that a helical filament with sufficiently large twist to bending rigidity and radius 
to pitch ratios, undergoes a sequence of stretching instabilities under the action of tensile force, even in the absence 
of any thermal fluctuations j^f. In the presence of thermal fluctuations or mechanical noise these instabilities are 
expected to lead to coexistence between different shapes of the deformed object [23|. 
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APPENDIX A: EULER ANGLES 



Since we are interested in three dimensions, we work in terms of the Euler angles 9, <p and ip. For convenience we 
will orient the force along the a;-axis, i.e. along the direction 8 — n/2, <fi = 0. The energy due to the force is then 

E F = -Fx = -F I dst% = -F I dssm(0(s))co$((p(s)) (Al) 
Jo Jo 

The Euler angles satisfy the differential equations pl|: 

d0 . , 

— = u>i sin ip + u>2 cos ip 

as 

sin 9 = —u>i cos ip + uj2 sin ip 

as 

— sin 9 — (ci>i cos ip — w% sin <p) cos 9 + lu^ sin 9 (A2) 
as 

We need to expand to quadratic order about the s = values 9q and </>o, and to linear order about ipo , (since the 
energy does not involve ip directly) giving 



d9 

— = K- X(lp - lp ), 

dep X + K(ip — ip ) cos 9q 



ds sin 8 sin 2 6>o 

dip cos #o 



X(9-9 ), 



ds sin #o 
where k and A are the rotated curvatures 



X + cj 3 , (A3) 



k = u>i sin ipo + u>2 cos ipo, 

X = —u>i cos -00 + ^2 sin ipo , ( A4) 



Thus, to second order in the curvatures 



9(s) =0 + [' ds' K (s') + ( f ds'X(s')) f ds'X(s') [ S ds"^(s"), 

Jo 2 sin O \Jo J Jo Jo 

<K*) = 4>o + -V T ds ' X ^ - ( f ds'K(s')) ( [ S ds'X(s') 

sin 9 J sin' 9 \Jo J \ Jo 



1 

sin 6*o Jo 



ds'K(s') / ds"iu 3 {s"). (A5) 
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APPENDIX B: FORCE-EXTENSION CURVE 



We will first compute the force- extensio n curve, via the partition function. We do this in two stages, depending on 
the value of the parameter \ = \JFL 2 /4a. 



Here we write 



where S is given by 



1. Small x (x *S 1) arbitrary fL) 



E F = —FL sin 9 cos (f>o + FS 



(Bl) 



S = — J ds cos 8 cos <t> J k — sin </>o y A — i sin # cos (j) (^j n 

— cos #o cos 4>o / ds'A(s') / LU3 — sin 0o / oVk(s') 

Jo Jo Jo Jo 

We expand e - FS / kT to quadratic order, yielding 



(B2) 



= y rf^o 



T>kT>\ exp 



f 2 ' 

1 - fS + —5 2 



(k 2 + A 2 ) + fL sin 9 cos </j 



(B3) 



where J dQ = sin 9 d9 d(j) Q . The functional integrals are all Gaussian and can be easily done by going to 

Fourier space. For example, 



so that 



L 

\Jo 



, , ^—v ,27ms, . ,27ms. 

K{S) = KQ + 2_^K n C0S(— — ) + }_^ K n Sm (— — ) 

n—1 n—1 



\ 2 1? ~ ( K c ) 2 + 3( K S ) 2 
oVk I ds = — Kq + i 3 2J n L,~ hr-^— + mixed terms 



(B4) 



3 " u ' ~ ^ 2(2nn) 2 

n=l K ' 

where the mixed terms will vanish latter upon functional integration and are omitted in the following. Defining 

and carrying out the functional integration by integrating over all the Fourier components 

00 00 



(B5) 



(B6) 



Vk = dn rj dn\ dn s j7 
i=i j=i 



we get 



where N — JT> k is the (normalization) integral over the exponential factor only. Similarly, 



(B7) 



(B8) 



J T>qk J (^J n^j ds = J Vqk 



KqL 2 



<L 2 



2 27m 

n=l 



= TV 



3a' 
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The functional integration with respect to A is carried out in the same manner. Putting this all together, we get 



fL 2 

1 — — — sin 9 cos < 



Z = N 2 j dn e fLsin8oCOs4 ' 
The integrals over solid angle are as follows: 

dn e fLs[n6 ° cos0 ° 
dQ a e fLsin00 cos *> sin6> cos O 



f 2 L 3 

(cos 2 9 cos 2 4> + sin 2 O ) 



6a 



sinh fL 



fL 

cosh fL sinh fL 



fL (fL) 2 
J tflo^-* ~*> cos 2 tf cos 2 O = 1 + jff/ £ - 

dO e /Lsinf, ° cos0 ° sin 2 0o 



Collecting all terms we finally get 



Z=AT' 



( sinh fL L 



6a 



cosh fL — 1 
(fL) 2 

sinh /L 



cosh /L 



(B9) 



(BIO) 



(Bll) 



From this, we can differentiate to get the average projection of the end-to-end vector on the d(\nZ)/df 



x = 



_L 1 _2)f 

tanh/L / 3/ 



fL 



tanh fL sinh 2 fL 



(B12) 



At large fL (but still small \), this goes as L — L 2 /(6a) — 1/f. For small /L this projection goes linearly with / as 
x ps /L 2 /3(l — L/(3a)) and we conclude that in the limit a» L the effective "spring constant" is 3kT/L 2 . 



1. Large /L (arbitrary x) 

We start by expanding 9o, 4>o about the values 0q, 0q which minimize Ef- To lowest order, 

'A 



y [ ds f ds') 
L Jo Jo 



so that the expansion reads 
Ef 



FL-H 
2L 



L ps \ j fL ps ^ 

ds / ds'n + I ds ds'X 
Jo I \Jo Jo , 



+ 



FL 



(9 - 9*) 2 + (0 O - 03) 



^lo'iilo^^Io^) 
Let us now consider the functional integral over k, the integral over A being exactly the same. We have 

Z K = J V Ke -§fo<*™ 2 -iSo Hlods'^T+MSo dsJ*ds' K ] 2 

= /^CXp{-^ K 2 -^f:(«) 2 + «) 2 ) 
J \ n=l 



2 



«) 2 + , «) 2 



2k k 



2 ^ 

12 ^V2(27rn) 2 2(2im) 2 (2nnf 



(B13) 



(B14) 



(B15) 
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We can diagonalize this expression by shifting defining 



with 



ah + b n 



b n = 



fL 3 



(27m) 5 



yielding an effective modulus of the zero mode 



,fL v-- "<> 



12 



ax 
tanhx 



We can now calculate Z K 



n=l 



aL 
2a n 



= (ao/a)' 

= {a /a)~ 
2X 



-1/2 



n 1 

n=l v 
-1/2 X 

sinhx 

1/2 



sinh 2\_ 

where Z k> q is the k partition function at / = 0. The entire partition function is then 



Z=Af 2 



2fL sinh2x 



For small x> this yields 



Z^M 2 



2 X 2 



which agrees precisely with the large / expansion of Eq. (|B11|) . 
The average extension in this second region (fL 3> 1) is given by 



X 



1 



2f /tanh2 X 



(B16) 



(B17) 



(B18) 



(B19) 



(B20) 



(B21) 



(B22) 



which for small \ gives x ~ L — 1// — L 2 /(6a), in accord with the large fL limit of the first region (\ <C 1). At 
large /£, it agrees with the Marko-Siggia result, Eq. I|14p, A uniform approximation for the force-extension relation 
is then 



1 



X 



2 X 2 ( 1 



fL 



2f /tanh2x tanh/L 3/ \t&nh fL sinh /i 



- 1 



(B23) 



This result is confirmed by simulations, as shown in data for a = 48, Fig. [21 It even works well when the persistence 
length slightly shorter than the chain length, as can be seen in the data for a — 6 in the right half of Fig. 0] Since 
our analysis is limited to the small fluctuation (stiff rod and/or large force) regime, we expect it to break down for 
very short persistence lengths in the small force limit. This is clearly observed in the case of a — 0.75, shown in the 
left half of Fig. 
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APPENDIX C: MODULI 



We now turn to the computation of the effective moduli of the various modes. We first treat the zero modes, 
starting as above with the small \ limit. We proceed by calculating a modified partition function, Z with a modified 
modulus of the zero modes, a, replacing the bare modulus a. Then 



L o 2 d\nZ 

2^ + A ° ) = -^r 



Tracking the changes induced in the partition function calculation above, we find 



-N 2 



sinh fL 



fL 2 









-J-Y 








- 


12a J 





giving 



so that 



lnZ = lna-lna + 21nAA + ln^/^-/L 2 ( -L 

fL \12a 



cosh fL sinh fL 

~jz (fiy 

1 \ / cosh/L 1 



24a J Vsinh/L fL 



(CI) 



(C2) 



(C3) 



(C4) 



(C5) 



(C6) 



At large /, the effective zero mode modulus was calculated above, Eq. (JB18I1 . The small \ limit of this agrees with 
the large force limit of Eq. (|C6|) . as it must. A uniform approximation for ap reads 



dlnZ 



da 

This implies an effective zero mode modulus of 

a = a + 

At small and large forces, this reads 

a 



1 

a 



12a 2 



cosh fL 1 



sinh fL fL 



fL 2 ( 1 



1 



12a Vtanh/L fL 



a + / 2 L 3 /(36a) fL < 1 
a + /L 2 /(12a) fL > 1 



a 



-fL 



1 



12a Vsinh/L fL 



tanhx 



(C7) 



This result is confirmed by the simulations as seen in Fig. El 

Now we turn to the finite wave-number modes. The interesting thing to note here is that the cos modes renormalize 
differently than the sin modes, due to the coupling of the zero mode to the cos modes. At large force, the effective 
modulus of the sin modes is the trivial 



fL 2 



(27m) 2 



(C8) 



so that the long wavelength modes are the most effected by the force. The effective modulus of the cos modes follows 
from 



(«) 2 > = (^) + 5<0 = ^- 



1 



with 



L(«) 2 ) 



(C9) 



(CIO) 



This expression is fairly messy, so let us examine its various limits. For large enough mode number, the effective 
modulus is a + fL 2 /(2irn) 2 as for the sin modes. For intermediate forces, all modes are large enough for this to be 
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L=6.283, a=48 




0.5 1 1.5 2 

F 



FIG. 6: Comparison of simulations for the zero-mode modulus as a function of the force F/kT to our uniform approximation, 
Eq. lICTt . 

the case. For strong forces, however, the effective modulus is L^/Ja/A up to some mode number, at which point it 
crosses over to the previous result. 
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